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Abstract 

A theoretical study is made on He scattering from three basic classes of dis- 
ordered adlayers: (a) Translationally random adsorbates, (b) disordered com- 
pact islands and (c) fractal submonolayers. The implications of the results to 
experimental studies of He scattering from disordered surfaces are discussed, 
and a combined experimental-theoretical study is made for Ag submonolay- 
ers on Pt(lll). Some of the main theoretical findings are: (1) The scattering 
intensities from the three disorder classes differ significantly, and can be used 
to distinguish between them. (2) Structural aspects of the calculated intensi- 
ties from translationally random clusters were found to be strongly correlated 
with those of individual clusters. (3) For fractal islands, just as for all surfaces 
considered here, the off-specular intensity depends on the parameters of the 
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He/Ag interaction, and does not follow a universal power law as previously 
proposed in the literature. 

In the experimental-theoretical study of Ag on Pt(lll), we use experi- 
mental He scattering data from low-coverage (single adsorbate) systems to 
determine an empirical He/Ag-Pt potential of good quality. Then, we carry 
out He scattering calculations for high coverage and compare with exper- 
iments for these systems. The conclusion is that the actual experimental 
phase corresponds to small compact Ag clusters of narrow size distribution, 
with partial translational disorder. 
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I. INTRODUCTION 



The properties of thin metal or semiconductor films on solid substrates are of major 
theoretical, experimental and technological interest. The structure of these films influences 
their physical and chemical properties, and these may be of central importance, for exam- 
ple, in the fabrication of electronic devices. The shape of adlayers generally depends on the 
growth kinetics and the microscopic characteristics of the substrate. Thin films are typically 
produced by epitaxial growth processes, wherein vapor atoms or clusters are deposited on a 
substrate.0'iiffii'Si Once adsorbed, the adatoms diffuse on the substrate surface, and when 
reaching short mutual distance they can form stable nuclei, which subsequently grow to 
clusters or islands by attachment of further adatoms. Nucleation and growth are competing 
processes, which depend on diffusion of adatoms, stability of the adsorbed clusters, and 
surface diffusion of these clusters. Therefore, thin film structures contain a wealth of in- 
formation about microscopic growth processes, and adatom/substrate and adatom/adatom 
interactions. 

Experimental studies have shown that in addition to ordered adlayers, various disordered 
structures may form. For example on adsorption of Ag on Pt(lll) surfaces, at low substrate 
temperatures and low deposition rates, isolated small clusters may form.iili'i Elevation of 
the surface temperature or lowering of the deposition rate can lead to the growth of fractal 
islands.ii At even higher temperatures randomly distributed compact islands may form.iili'i 
Therefore, thin films offer an exceptional opportunity for studying two dimensional disor- 
dered systems. By investigating the relation between the structure of thin films and their 
growth processes, valuable information can be obtained about the processes leading to the 
emergence of different kinds of disorder. This, however, requires techniques for determining 
the disordered adlayer structure. 

Surface morphologies can be probed by scanning tunneling microscopy (STM) or diffrac- 
tion techniques, such as thermal energy atom scattering (TEAS), low-energy or reflective 
high energy electron diffraction (LEED or RHEED). Diffraction methods have the property 
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that in addition to being sensitive to local surface features, they resolve the global sur- 
face structure. Accordingly, they provide a better average of the over-all topography than 
direct imaging techniques, such as STM. Moreover, diffraction readily allows for monitor- 
ing adlayer growth in situ and at various temperatures and therefore can easily follow its 
temporal evolution. For studying surface structure, a technique of great power is Helium 
atom scattering. It is non-destructive, and He only probes the outer layer. Also, thermal 
He atom scattering is dominantly elastic. Moreover, since its wavelength is of the order 
of surface unit cells, it is very sensitive to the local adlayer structure. He atom scattering 
has successfully been employed in the study of ordered surfaces and surfaces with isolated 
defects. However, relatively little is known as yet on the manifestations of different kinds of 
disorder in He scattering patterns. Two of the few experimentally and theoretically stud- 
ied disordered surfaces are those of substitutionally disordered mixed Xe+Kr monolayers 
on Pt(lll),i and translationally disordered small Pt-clusters on Pt(lll).0 It was found 
that attenuation of the specular peak in He scattering from the surface due to the presence 
of adsorbates or other defects contains a wealth of information about adatom/adatom in- 
teractions, the clustering of adatoms, and 2D vs. 3D epitaxial growth.0 Also, non-Bragg 
maxima, e.g., Fraunhofer or rainbow maxima, appear for such systems. These maxima con- 
tain information on the microstructure of islands and are sensitive to percolation transitions. 
Non-specular Bragg peaks appear as well, which contain information on the corrugation of 
the He/adlayer interaction.! The study of randomly adsorbed rare-gas overlayers has the 
advantage that the interactions between the rare gas atoms and He are well known. In 
contrast little is known about the interactions in a He/ (metal adsorbate)/ (metal substrate) 
system. So far, this has complicated the interpretation of features in the angular scattering 
patterns, and has made identifying different kinds of disorder difficult. 

In this paper, for the first time a comparative He scattering study is undertaken among 
three major classes of disorder: Translationally random adatoms, translationally random 
compact islands, and fractal/dendritic submonolayers. These classes are representative 
of a large proportion of the experimentally reported adlayer structures at submonolayer 
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coverage.0 It is thus of fundamental interest to understand the differences in their respec- 
tive He scattering intensity distributions, and to investigate whether the type of disorder 
can be identified by a He scattering experiment. In order to carry out this program rehab ly, 
we determined the interaction potential between He and a Pt(lll) surface with an adsorbed 
Ag atom. From this we have constructed what we hope is a realistic potential function 
for He interaction with a disordered metal submonolayer. The choice of the Ag/Pt(lll) 
system was motivated by the fact that the He/Pt(lll) potential is rather smooth parallel 
to the surface. Therefore the structures in the scattering patterns are mainly due to the Ag 
adsorbates. Moreover, the He/Ft (111) belongs to one of the few atom/surface systems for 
which a fairly reliable empirical potential is available. 

The outline of the paper is as follows: Section H describes the experimental and theo- 
retical methods employed, and states the assumptions made in the theoretical part of this 
work. Section |IT| discusses the scattering from the different disorder classes. In section [HI C 



we analyze the experimental scattering data and suggest a morphological identification of 
the disordered phase present on the surface. Concluding remarks are presented in section 

m 

II. SYSTEMS AND METHODS 

The adlayer structures studied here theoretically are: 

1. Full Ag monolayer. 

2. Single adsorbed Ag atom. 

3. Isolated small compact cluster (henceforth SCC). 

4. Ag adatoms randomly adsorbed on a lattice. 

5. SCCs of Ag randomly adsorbed on a lattice. 

6. Ag adatoms completely randomly and continuously distributed (Ag adatoms can 
be arbitrarily close). 



7. Large compact Ag islands. 

8. Fractal islands of Ag atoms. 

(|T]) and (0) are important as a reference, and for developing potentials. All other cases 
were chosen as representing important types of disorder. In all cases, the adsorbates are 
on a Pt(lll) surface. In the case of randomly adsorbed Ag adatoms one would expect to 
observe system and not the idealized system (P). However, calculations for the latter 
structure leads to additional insight into angular scattering patterns. 

A. Experimental Methods 

The experiments were performed in a high resolution UHV helium scattering apparatus 
with a nominal base pressure < 10^^° mbar. The system is equipped with a commercial 
Knudsen cell by means of which high purity (99.999%) silver can be evaporated onto the 
Pt(lll) sample at a rate between 10^^ and 10^^ monolayers per second. The Pt(lll) sample 
was cleaned by repeated cycles of heating in an oxygen atmosphere (10~^ mbar) at 750K, 
sputtering with 1 keV Ar-ions followed by short annealing at 1250K. The sample is mounted 
on a manipulator which allows the crystal to be rotated around its three principal axes. In 
addition, the total angle x = + ^/ of the incident and outgoing He-beam can be varied 
between 60° and 120° by rotating the detector. By means of a liquid helium cryostat the 
sample can be cooled to 20K and by concomitant heating by electron bombardment the 
surface temperature can be set and held constant at any value between 20K and 1250K. 

The experiments reported here (e.g., Fig.|l]) involve two different modes of operation of the 
He-scattering apparatus. In the first mode, the specularly scattered He-intensity is recorded 
during the deposition of silver at a constant rate. From the initial slope (extrapolated 
towards silver coverage — * 0) the scattering cross-section of the Ag-nuclei initially formed 
on the surface are obtained. From a detailed analysis of the shape of the deposition curves, 
we can infer that at the temperature of 38K at which the data in Fig.|^ were taken these 
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initial Ag nuclei are individual Ag atoms randomly distributed on the Pt(lll) surface. The 
initial slope of the He- intensity curve, therefore, directly yields the He cross section ^(A;^) 
(as defined by Eq.(|^) below), where kz = k cos(^j) is determined by the wavevector k of the 
incident He beam, related to its energy by E = h'^k'^ / {2m). 9i = Of (specular scattering) 
denotes the incident and outgoing angle of the He beam measured against the surface normal. 
To measure the dependence of the cross section S on k^, either the total scattering angle 
X = Oi + 6 f OT the energy of the incident He beam has to be varied. We have chosen the 
second option which, in practice, we achieve by varying the nozzle temperature from liquid 
nitrogen temperature to above room temperature at constant total scattering angle x = 90°. 
This temperature variation corresponds to a change of the energy E from 18 to 90 meV, 
i.e., from k^ = 4.151"^ to 9.281"^ 

The second type of experiment is the analysis of diffraction profiles as shown, e.g., in 
Fig.^]. Here, the He-intensity is recorded as a function of the parallel wavevector transfer 
AK. In the present case, the incident and outgoing beams lie in the same plane as the 
surface normal {in plane scattering) and AK = k{sin6f — sin^j). The azimuthal orientation 
of the scattering plane relative to the surface crystallographic directions can be varied by 
rotating the crystal around its surface normal. In Fig.^ the surface was oriented along the 
[112]-direction (i.e., the direction along which the first order Bragg peaks of the hexagonal 
substrate lattice are expected). 

Our scattering apparatus is equipped with a time-of-flight (TOF) spectrometer allowing 
to separate the elastic from the inelastically scattered He intensity by counting only those 
He-atoms within a narrow energy window of about O.SmeV centered around the incident 
beam energy. In this way, even the very small elastic signals away from the Bragg peaks can 
be discriminated. As shown in the next sections, these elastic features provide important 
information on the surface morphology. 

All the experiments reported here were conducted at low temperatures, T = 38K. As 
already mentioned, this ensures that small silver clusters are formed as a consequence of the 
limited mobility of the Ag adatoms on the Pt(lll) surface. Although similar in size and 
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distribution, these low-temperature structures should not be confused with those obtained 
after depositing or annealing Ag at surface temperature, T > 620K.I'II1 At these elevated 
temperatures, the Ag atoms are embedded into the topmost Pt(lll) layer forming small, 
stable clusters as a consequence of the surface strain.iS. Instead, the low temperature 
structures presented here are diffusion limited aggregates composed of Ag atoms adsorbed 
on top of the Pt(lll) surface. A detailed account on the difference of these two configurations 
will be given in a forthcoming paper. 



B. Theoretical Methods 

In modeling the substrate, we assumed a flat supporting surface. Indeed, the He/Pt(lll) 
equipotential surface is rather smooth, leading to very little off-specular scattered He.0 How- 
ever, due to the lattice misfit of about 4% between Ag and Ft, the assumption of a flat 
Pt(lll) surface is only a reasonable first approximation. In addition, in all of our calcula- 
tions we assume the Ag/Pt(lll) system to be rigid. Note that experimentally the inelastic 
contribution can be separated by a TOF spectrometer. As for elastic scattering, results of 
calculations using a rigid, non- vibrating surface system should be useful at least for studying 
the main qualitative features. Moreover, for diffraction scattering from crystalline surfaces, 
the effect of surface vibrations on the scattering intensities can be represented approxi- 
mately by a simple Debye- Waller factor. A similar description should be successful for 
the angular intensity distributions obtained in scattering from an adlayer. 



1. Kinetic Monte Carlo Simulations of Diffusion Limited Aggregates 

We outline here the methods used to simulate the compact and fractal Ag islands (sys- 
tems Figs.^l^. There is an extensive literature on the modeling of adsorption, diffusion 
and aggregation processes on homogeneous surfaces.^ To produce realistic arrangements of 
adatoms on the surface we adopted the Kinetic Monte Carlo (KMC) method.i2l The Ft (111) 
surface was simulated as a hexagonal lattice of 100x100 unit cells with periodic boundary 
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conditions. As for the lateral interaction between adsorbed particles, the assumptions made 
are critical for the shape and size of the final clusters formed. We assumed that the interac- 
tions between two adsorbates depend on the number of their nearest neighbors, tinn- Thus 
the energy of a particle in a given configuration is: 

E(n) = Eo + n^N e (1) 

where Eq = 5.2kcal/mole is the activation energy for diffusion at zero coverage, and e = 
5.0kcal/mole is the nearest neighbor interaction energy. No values were available for the e 
parameters for Ag on Pt(lll). Therefore, we adopted the values given by Rosenfeld et al.i 
for Pt adsorbed on Pt(lll). These values, while not more than semiquantitatively valid, 
should at least allow for the formation of plausible aggregates. 

The diffusion of adatoms was modeled by random walks over nearest neighbor sites. The 
hopping rate, i.e., the transition probability per unit time, of a particle from site i, where it 
has Hi nearest neighbors, to a neighboring site with Uf nearest neighbors, is taken as 

Ui^f = zye-^^"')/^'^^) = uoe-^''/^''^^ (2) 

where uq = u exp[— _E'o/(^^)] is the hopping frequency of an "isolated" particle; tq = 1/ujq 
is the average time interval between successive moves on the bare (zero coverage) surface. 
Since there is an uncertainty in our data regarding the values of the activation energy Eq 
and the interaction energy e, it must be commented that there is a resulting uncertainty 
in the temperature scale: Changing T at constant E{n) is equivalent to changing E{n) at 
constant T. Note further that according to Eq.(^, the transition rate depends only on 
the initial state. A possible dynamical interpretation of this model, originally suggested by 
Uebing and Gomer,ill is that the rate of the transition from state i to state /, is governed 
by the energy difference between the initial state and the transition state. The same is of 
course true also for the reverse process (with the final state replacing the initial state), so 
that detailed balance is obeyed. We note that a situation where a rate depends on the gap 
between any given initial state and the transition state holds widely for many activated 



rate processes. Hopping diffusion involves overcoming a barrier between the initial and final 
configurations and therefore the model seems most reasonable here. The KMC simulations 
were performed using the time- dependent Monte-Carlo scheme.^ In this scheme, instead 
of randomly choosing particles and accepting or rejecting moves according to the given 
transition probabilities (as in the "traditional" MC simulations), one performs a move in 
any attempt, and propagates the time accordingly. More explicitly, we first calculate the 
average transition rate (r) = 1/ (uJi^f) out of state i, then randomly sample a given i ^ f 
move with probability uJi-^f / J2 ^i^f and, finally, perform this move and record the time 
elapsed as At = 1 / (r) . Further details about this procedure can be found elsewhere.E10 

The simulation starts with two particles diffusing on the surface until a third particle 
adsorbs. This period is typically 0.05 sec (for simulating an adsorption rate of lML/500sec), 
corresponding to about 10^ MC steps at room temperature. Then the three particles diffuse, 
a fourth one adsorbs, and so on until the desired coverage is obtained. The configurations 
attained are not at equilibrium, but reach a steady state after some time duration, when 
the adsorption is terminated. The true equilibrium state corresponds most probably to a 
segregated phase, with one large compact island formed by all the adsorbates. Both the 
compact and the fractal clusters of Figs.^J^ were generated according to this procedure, 
differing only in the respective temperatures of 500K and 200K. 



2. He Scattering Calculations 

The scattering intensities were calculated using the Sudden Approximation 
which has proved very useful in studies of atom scattering from defectfl 
and from substitutionally disordered rare-gas monolayers.! On the basis of the experience 
gained with the SA, including tests against numerically exact calculations for several model 
systems,^ we estimate that at least the main predictions of the SA calculations should be 
reliable for the systems studied here. 

The SA takes the following form for the systems studied below. Consider a particular 
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(static) configuration r = {fi, ...,r/v} of a disordered adlayer system consisting of atoms 
with fi denoting the position of atom i. The angular intensity distribution for He scattered 
from this adlayer involves an average over all configurations pertinent to the disorder:iiii 



(3) 



In this expression AK = K' — K is the wavevector transfer of the He parallel to the surface, 
where {K, kz) is the incident wavevector, and {K', —kz) the final wavevector of the scattered 
He atom; R = {x, y) denotes the coordinates of the He atom in the surface plane; A is 
the area of the surface over which the integration in Eq.@ is performed; (...) denotes the 
average over all configurations of the quantity in parenthesis; rif{R) denotes the phase shift 
for He scattering for an adlayer having a configuration r. The phase shift is given in the 
WKB approximation by@ 

7^^{R)= I dz ( kl-2mVriR,z)/h^ -kz) -kzi{R), (4) 

where m is the mass of the He atom and z denotes the distance of the He atom from the 
surface plane; V5?(-R, z) is the interaction potential between the He atom at position {R,z). 
^ (R) in Eq. (|^) is the classical turning point for the He atom when colliding with the surface 
at the lateral position R, i.e., ^(-R) is the z value for which 

{hkzf - 2mV^{R, z) = (5) 



3. The Interaction Potential 

To determine the interaction of He atoms with the Pt(lll) surface with a single adsorbed 
Ag atom, we calculated the attenuation of the specular peak with incident He energy. This 
attenuation is closely related to the cross-section of the adatom. The relation between the 
cross-section S and the specular intensity I for scattering from random isolated adatoms 
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where Iq represents the specular scattering intensity from the clean Pt(lll) surface, n is the 
number of sites per unit cell, and O the adatom coverage. The parameter S in Eq.(P) can 
be interpreted as the cross-section of a single adatom. Indeed, Eq. is only valid when the 
cross-section of different atoms can be viewed as non-overlapping. S is very sensitive to the 
interaction potential. Therefore, we tried to find a potential that is capable of reproducing 
the experimentally observed changes in the specular scattering. We assumed the following 
general form for the interaction potential between He and the Pt(lll) surface with an 
adsorbed Ag atom: 

^He/Ag/Pt(lll) = VHe/Pt(lll) + ^He/Ag, (7) 

where VHe/Pt(iii) is the interaction potential of He with the clean Pt(lll) surface. A He-Pt 
Morse potential, 

V{z) = De (^-^'"^ - 2e-" (^-^'")) (8) 

for this interaction^ is available in the literature, with De = 7.86meV= 2.89 x 10~^a.u. and 
a = 0.98A"^ = O.52a.u.0 Vne/Ag is the interaction between He and an isolated adsorbed 
Ag atom. We modeled this interaction with a Lennard- Jones potential, and changed its 
parameters until the potential could reproduce the dependence of the experimental cross- 
section on the incidence energy of the He atom. This procedure also yielded the value of 
Zm = 11.46a.u. in the Morse potential. The high sensitivity of this fit to the parameter 
values is the source of our confidence in the semiquantitative reliability of the potential. 

For the interaction between He and multiple adsorbed Ag atoms, we assumed pairwise 
interactions for the He atom with the adsorbed Ag atoms, i.e., Vne/Ag in Eq.(|^) is replaced by 
^iVue/Agifi) where rj is the distance of the He from the i*^ Ag atom. The pairwise additive 
potential for the interaction with a collection of Ag atoms is expected in general to be no 
better than a first approximation. We estimate, however, that it is sufficient in accuracy for 
our purpose here. 
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III. RESULTS AND DISCUSSION 



A. He Interaction Potential with Pt(lll) and adsorbed Ag 

Specular scattering intensities were measured as a function of He incidence wavenumber 
for a single adsorbed Ag atom. Using these values in Eq.(^ provided us with the exper- 
imentally determined cross-sections, shown in Fig.|l|. Experimental constraints detailed in 
Sec. [ll A| limited the range of incident He wavenumber values. Subsequently, we attempted 
to fit the experimental cross-section data by finding an optimal set of C12 and Cq values for 
the He/(adsorbed Ag atom) Lennard- Jones 6-12 potential 

Vue/Ag - ^ - 

with which the specular intensities were calculated from Eq.(P). By adopting this trial 
and error approach, we found that the experimental values are best reproduced with the 
following Lennard- Jones parameters: C12 = 6.33 x lO^a.u. and Cq = 29.0a.u. This leads to 
reg = (2Ci2/C6)^/^ = 4.6A and Vir^,) ^ 0.9meV. A comparison between the experimental 
and theoretical values is presented in Fig.|I[ We found the fit to be very sensitive to these 
parameters, suggesting that our empirically determined potential models the long range 
attractive and the short range repulsive interactions well, and its accuracy is mainly limited 
by the experimental uncertainty. 

B. Angular Scattering Patterns from Disordered Overlayers 

The results will be outlined by examining the angular distributions of each of the disorder 
classes, and comparing among them. In each case where random ensembles are discussed, 
we averaged over 30 configurations. 
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1. Scattering from a Full Monolayer (system ^) 



To evaluate the influence of surface disorder on He scattering, we flrst calculated for 
reference the angular intensity pattern for He scattered from a completely ordered monolayer 
of Ag on Pt(lll), i.e., for a silver coverage = 1. (The issue of incommensurability 
does not arise here because we take the Pt(lll) surface to be flat - Eq.(^)). The incident 
perpendicular He momentum is 3bohr^^ = 5.Q7A^^ (For the other calculations presented 
below, we assumed the same incident momentum.) The result, presented in Fig.|^, strong 
flrst order Bragg peaks along with the specular peak. The very weak background intensity is 
due to flnite grid effects. Note that our calculations, based on a pairwise additive potential, 
probably give the effect of a too corrugated surface (hence diffraction stronger than the real 
one) since the smoothing role of the conduction electrons is not included in the interaction. 

2. Scattering from Single Clusters: Adatom and Small Compact Cluster (systems ^,^) 

The simplest deviation from perfect order is the introduction of a single defect. Such 
defects, namely a single Ag adatom and SCC are the subject of this section. In subsequent 
sections gradually more disordered systems will be considered. 

The dashed lines in Figs.|^,^ present the scattering results from a single adatom and 
an SCC (heptamer) respectively. Despite their quantitatively different appearance, the 
scattering patterns share some important qualitative features: 

a. Specular Peak: The sharp specular peak is due to He scattered from uncovered 
Pt(lll) surface areas, where the He/Pt(lll) potential is not influenced by the defect. The 
flatness of this potential across the unperturbed surface results in almost perfect mirror 
scattering. 

b. Broadening of the Specular Peak: Immediately below the specular "spike" a shoulder 
appears, which is clearly wider in the single adatom case. In other words, the broadening of 
the specular peak is inversely related to the cluster size. 
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c. Off-Specular Peaks: A series of rather broad, off-specular diffraction features follows 
the specular broadening. The interesting feature is their spacing, which is larger for the 
single adatom. Again, this spacing in reciprocal space is inversely related to the cluster size. 
For a detailed quantitative analysis of this issue see Ref.||TO[. Here it will suffice to mention 
that the peak structure, as well as the specular broadening, can be satisfactorily explained in 
terms of scattering from a hard hemispherical object of radius d. The resulting Fraunhofer 
diffraction intensities are given byS 

'1 + COS0) Ji{kd sm(f:^ 



/(0)oc 



(10) 



sm(p 

where (p is the scattering angle, Ji denotes the Bessel function of first order, and k is the 
incident wavenumber. The parameter d can be interpreted as the radius of the scattering 
cross-section. Indeed, we found using Eg. (p!OD , that for scattering from an isolated heptamer 
(Fig.^), the Fraunhofer interferences are those characteristic for an average island scattering 
cross-section of three nearest-neighbor Ag atoms. 

Another mechanism at work in the scattering from defects is rainbow scattering, in 
which the positions of the peaks are determined by the inflexion points of the equipotential 
surface.ii'Ei For an extensive treatment of this issue in the heptamer case, see Ref . ||10[ . 



3. Scattering from Adatoms and SCCs Randomly Adsorbed on a Lattice (systems ^,^) 

A natural transition to more complex disorder is obtained by considering a collection of 
clusters, randomly adsorbed on the hexagonal lattice. A system like this can be formed in 
realistic conditions when gas phase Ag atoms are deposited on a cold Pt(lll) surface, such 
that the Ag mobility is very small.i'i0i Scattering from such translationally random systems 
is the subject of the present section. 

The results for a system of single adatoms and SCCs at 15% coverage are presented as 
the solid lines in Figs.|,|[ Let us enumerate the central features: 

a. Specular and Bragg Peaks: The sharp specular peak is as usual an indication that 
significant portions of the Pt(lll) surface remain flat. In addition one observes strong 
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first and second order Bragg peaks at multiples of 2Ti/a (a = 2.77 A is the Pt(lll) lattice 
constant). What is their origin? It is not the underlying flat Pt(lll) surface, as proved 
by the absence of Bragg peaks in the case of the individual defect systems of Figs.^,^. 
Nor is it the small (111) plateau on top of the SCCs (see inset of Fig.^), since the Bragg 
peaks are observed also for the adatom system, which has no such structure. Thus these 
peaks can only be the result of the strong corrugation induced by the presence of the Ag 
adatoms and clusters on the lattice. This must be so in spite of their random positions on 
this lattice. Since the Bragg peaks are robust and their intensities are comparatively high 
(~ 10"'^ relative to specular), they are easily measurable. In addition, they yield to a simple 
theoretical description, and thus are attractive candidates for the calibration of a simple 
model potential from experimental data. 

b. Similarity to Single Cluster and Effect of Disorder: Perhaps the most striking result 
from the present calculations is the close (qualitative) resemblance between the intensity 
distributions from the single clusters and their random counterparts. This can be seen 
very clearly in Figs.|],^, where the distributions are slightly offset for clarity. As is especially 
conspicuous in the SCC case (Fig.^), there is a one-to-one correspondence between essentially 
every (broad) off-specular peak in the single and random case, respectively. The effect of the 
randomness seems to be limited to a modification of the single-cluster intensity distribution 
by the addition of noise. In the previous section we identified the principal cause for the 
peak structure in the scattering from a single cluster to be Fraunhofer scattering. This, 
then, must also be the dominant mechanism in the scattering from the random systems. 
The Fraunhofer mechanism is sensitive to local surface details, through the cross-section of 
the individual cluster [Eq.(|lOD]. Since this is a feature of the adlayer which is not expected 
to change significantly in the transition to an ensemble, the similarity agrees with intuition. 
But why does the randomness have such little effect? This can be understood qualitatively 
as follows: For a regular super-lattice of islands, one would expect a set of Bragg peaks at 
multiples of the inverse super-lattice constant. For a random set of islands, as in the present 
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case, such long-range translational order is not present, and the only order that remains in 
the system is the structure of the individual islands and the discreteness of the underlying 
lattice giving rise to the Bragg peaks at the Pt positions. Hence this is the only coherent 
contribution (though the Fraunhofer or rainbow mechanisms) to the intensity spectrum. 
The effect of the disorder is then reduced to the addition of noise. In the next section we 
investigate what happens when even the pinning due to the underlying substrate lattice is 
lost. 

4- Scattering from a Fully Random Adatom System (system ^)) 

In the previous section we considered a random adatom system on a lattice. Here we 
consider the same system of randomly located adatoms, with the lattice constraint removed. 
Thus, each Ag adatom is located completely randomly, as if the Pt(lll) surface were per- 
fectly smooth. Although this system is rather artificial, it removes the distance length scale 
between the atoms, which the lattice implies. Consequently in this model, adatoms can be 
"merged" into each other into statistically allowed structures. This non-lattice system is 
amenable to analytical modeling, as discussed below. 

The scattering results from this system are presented in Fig.H, for a system of 15% 
coverage Ag on Pt(lll). The scattering intensity is very similar to that obtained for the 
lattice system of Fig.^ The broadening of the specular peak is essentially identical, and the 
decay of the intensity beyond the broadening is the same apart from the oscillations absent 
in Fig.^. These oscillations, then, can most likely be identified as the result of Fraunhofer 
interference, whose effect is suppressed due to the loss of adatom identity in the present 
case. 

Another sharp feature evidenced in Fig.^ is the presence of two distinctly different slopes 
dominating the intensity distribution. The slope discontinuity occurs at about /\K* = 3A~^, 
and is an indication that there exist two important structural regimes on the surface: The 
individual cluster, and the long-range translational randomness. Why SA^^? The answer 
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can be found in the single adatom scattering intensity of Fig.^: 3A~^ is exactly the point 
where the specular broadening reaches a minimum, and the Fraunhofer interference pattern 
takes over (the actual value of course depends on the cluster radius). One would expect that 
the small AK (< AK*) regime corresponds to the long range surface-structure, and vice 
versa. To check this, we repeated the calculation with a "harder" He/Ag potential (C12 = 
3.0 X 10 ''a. u., Cq = T.Oa.u.), which enhances the effect due to the individual cluster structure. 
The result is presented in the right inset in Fig.^. The slope discontinuity at AK* is even 
more pronounced here, and oscillations reminiscent of the on-lattice system are observed as 
well for AK > AK*. The reason is that this harder potential effectively increases the local 
surface corrugation; consequently the surface is rougher and leads to increased high-angle 
scattering. Thus changes in the He/Ag potential strongly affect the intensity distribution 
for large AK values, which therefore contain important information on the He/adsorbate 
potential. Furthermore, by studying the large AK regime, a He scattering experiment can be 
used as a tool to measure adatom electron densities, which are related to the He/adsorbate 
potential.ii It is clear from the above analysis that the short-range structure, manifested by 
the local corrugation (and controlled by the hardness of the potential), appears in the large 
AK regime. However, the small AK behavior is identical even for the hard-potential case. 
This proves that this regime is determined exclusively by the long-range structural features 
of the surface, as indeed expected. 

The question arises whether we can isolate the long-range effect of the surface disorder 
{AK < AK*) from that of the interaction with a single defect. In order to completely isolate 
the effect of the translational randomness of the clusters, we introduce a highly simplified 
two-state model (TSM), in which every adatom is represented by a cylinder of height h 
and a diameter d, and the adatoms are completely randomly located. Here h is the height 
of an adatom above the Pt(lll) surface and d is the diameter of the adatom scattering 
cross-section. The adatoms are assumed to be fully mutually penetrable. The resulting 
system may be considered as a two-phase random medium, in the sense that the underlying 
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surface and the incomplete layer of tops of islands are two separate phases from the point 
of view of the incident He beam. This approximation is reasonable for systems of metal 
atoms adsorbed on metal surfaces, because the corrugation of metal islands is extremely 
weak. By representing the adatoms as structureless cylinders, we capture only the global 
features of the system, neglecting its local structure, and can thus hope to reproduce only 
those features of the spectrum which result from the underlying random scattering system. 

An outline of the solution of this TSM for the intensities is given next; for additional 
details see the appendix. For scattering from a hard island on a hard surface, the general 
expression for the scattering intensity [Eq.(|^)] becomes: 

j(Ai?) = l^ldR, I dR2 e^A^-(^i-^2)(e2^'=^^(^^'^^)), (11) 

where 

Z(i?i,4)=^(i?2)-^(i?i) (12) 
and ^(-R) is the classical turning point at R for a hard surface, and 

^^2ik.ziR,A)^ = J dze^"'^'fz{z;Ri,R2). (13) 

Here fz{z; Ri, R2) is the probability density of observing a height difference z between two 
points located at Ri and R2 on the surface. The problem of obtaining I{Q) therefore reduces 
to the evaluation of fz{z', Ri, R2)- In the case of realistic shape functions for the island, the 
evaluation of fz is a very hard problem, because it requires to take into account the 3D shapes 
of all possible cluster types. The TSM renders the problem of finding fz tractable, at the 
expense of loosing certain details in the long-range part of the real intensity distribution. 
Rewriting Eq.(|T^) for the case of a stepped surface, we obtain: 

(^^ik,z(A,R2)^ = Y^e^ik,'^Vi{zf,Ri,R2) (14) 
i 

where the sum is over all possible height differences, and Vi^Zj] Ri, R2) is the probability 
to observe a given height difference between the points Ri and i?2- In the present case 
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Zj = 0,±h. The problem now reduces to evaluating Pt^zj; Ri, R2). Our evaluation for 
randomly adsorbed islands is presented in the appendix, and shows that for this case the 
scattering intensity Eq.(|l^) becomes: 

I{AK) = Clip d"^, 2k, h,A)6{AK) + C2ipd^,2k,h,d^/A) 

gia,(3)= t UaS)e^^^'°'~"^-'^'^^)sdS, (15) 
Jo 

where Ji is the modified Bessel function of order 1, and p is the adsorbate density. The 
factors Ci, C2 are independent of AK. The important point to notice about Eq.(|l^) is that 
apart from a specular term, it is composed of a term due to a single cylindrical adatom (the 
J I term - compare to Eg . ([T0|) ) , and a new term ((?), which incorporates the effect of the 
surface disorder. As can be seen in the left inset in Fig.^, the effect of the disorder term 
is very significant: It completely dominates the small AK regime compared to the single 
adsorbate term. In particular the slope in this AK range is controlled almost exclusively by 
the disorder term (in contrast to the width, which was previously shown to be determined 
by the radius of the islands). This explains the universality of this slope in all the random 
models considered above. 

Next, Eq.(|I^) was applied to study the quantitative aspects of the small AK region. 
The height h has mainly the effect of an overall intensity factor. Since one cannot expect 
the crude TSM to correctly predict the intensities, we shifted the intensity curves in the 
(Ai^', log(/)) plane, and focused on the effect of changing the cylinder diameter d. At 
the incidence wavenumber of = 3.0bohr~^ = 5.7A~^, Fig.[l| predicts a cross-section of 
~ 140A^, i.e., a diameter of ~ 5a. In Fig.|] we show the TSM prediction for d = 5a, 10a. 
The effect of increasing the diameter is to flatten the intensity. The d = 5a case clearly 
underestimates the intensity, whereas for d = 10a there is good quantitative agreement in the 
small AK (< AK*) regime. Thus, not surprisingly, the TSM does not offer good quantitative 
agreement for realistic parameter values. However, the agreement improves when the TSM 
is compared to the "hard-potential" case (see inset in Fig.^), which is consistent with its 
formulation hard-wall model. 

20 



In conclusion, the relative success of the TSM in representing the small AK regime 
(as well as its failure to agree with the large AK region!), demonstrates convincingly that 
the properties of the long-range structure of the disordered surface can be found from the 
quantitative analysis of this regime. 



5. Scattering from Large Compact Islands ( system ^ ) 
In this and the next section we discuss the scattering from yet more complex, stochastic 



systems: those created by the Kinetic Monte Carlo algorithm described in Sec. |IIB 1| . A 
typical configuration of the large, quasi-compact islands formed at 500K is shown in the 
inset in Fig.^. Note that periodic boundary conditions are employed. It is not our purpose 
here to obtain a quantitative description of the structure of these islands from the scattering 
intensities. Instead, we will limit ourselves to a discussion of some of the main features 
of the intensity distribution, and contrasting them with those characteristic of the systems 
studied above. 

a. Broadening of the Specular Peak: Among all the systems discussed so far, the spec- 
ular broadening observed in Fig.^ is the smallest. This corresponds, as expected, to the 
large area of the individual islands. The observed width is ~ 0.3A~^, in agreement with 
an average island diameter of 13-17a obtained in our KMC simulations. Thus, as indeed 
obvious, the specular broadening can be used experimentally to find the average diameter 
of randomly adsorbed islands. 

b. Broadening of the Bragg Peaks: In contrast to the case of translationally random 
adatoms or SCCs, there is little off-specular structure in Fig.||. The exception are the 
Bragg peaks, whose existence in the present case is not surprising, considering the large 
island surface area, which has the underlying Pt(lll) lattice corrugation. However, the 
new feature is the broadening of the Bragg peaks. The broadening is roughly half that 
of the specular peak. This phenomenon is reminiscent of Fraunhofer interference of light 
transmitted through a grating. Indeed, roughly speaking, the uncovered surface area between 
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the islands can be thought of as a sht pattern. The kinetics governing the formation of the 
islands leads to depletion zones, or an effective repulsion between them (since they cannot 
overlap), thus inducing a degree of regularity in this pattern, which can yield the Fraunhofer 
broadening. Note that this is different from the Fraunhofer scattering by a hemispherical 
object described by Eg. (]T(]|) , which is responsible for the specular peak broadening and the 
"humps" at AK = ±1.5, 3.0A~^. The irregularity in "slit" sizes and orientations, as well 
as to a lesser extent their positions, both due to the inhomogeneity of island sizes and 
separation, is responsible for the added noise. 

6. Scattering from Fractal Islands ( system ^ ) 

In this section we discuss the scattering results from the last model system considered 
in this work: the dendritic/fractal structures formed by the KMC simulations at 200K. The 
interest in such structures hardly needs recapitulation. The averaged scattering intensities, 
as well as a typical configuration, are shown in Fig.|^. This configuration is rather reminiscent 
of STM results on Pt/Pt(lll)0ii and Ag/Pt(lll)S obtained under similar conditions. The 
fact that these systems are indeed fractal over a number of orders of magnitude of resolution- 
observation, is discussed in detail elsewhercil 

Interestingly, in spite of the non-negligible surface area of the clusters, there is apart from 
weak first order Bragg peaks almost no trace of interference in Fig.|^. This is in agreement 
with the characterization of a fractal as an object with no typical length scale between 
physical lower and upper cutoffs. Indeed, our "fractals" are composed of a large variety of 
clusters of different sizes (dimers, trimers and larger clusters), which generally are somehow 
connected. Each of these clusters contributes Fraunhofer peaks at its typical length scale. 
The net effect is to smear out these peaks into a relatively smooth decline of the intensity. 
However, since the fractal clusters are all finite in size, their average radius is reflected in 
the scattering distribution, as the by now familiar specular broadening. Apart from this, the 
absence of a typical length scale (above the lower cutoff a) results in that there is essentially 
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no structure left in the intensity distribution. 

As for the intensities at beyond the broadening, we observe that the dechne of the 
angular intensity with Ai^ for fractal islands is very similar to that for randomly adsorbed 
adatoms (Fig.|^), and quite similar to the average decline for SCCs and large compact is- 
lands (Figs.^J^). Furthermore, as mentioned earlier, this dependence is determined by the 
He/Ag/Pt(lll) interaction potential. Nevertheless, we attempted to fit the intensities both 
before and beyond the broadening (|Aif| > 2Ar^) by a potential- independent power-law, 
following Ref . |36| : 



/(AK) oc Air^-^ (16) 

Here D is the fractal dimension. The right inset of Fig.|^ shows the result of this fit for both 
positive and negative /S.K (the latter were reflected), after elimination of the specular and 
Bragg peaks. The result is a slope of -1.84 ± 0.03 for \/\K\ < 2l-\ and -1.62 ± 0.06 for 
|Ai^'| > 2A~^ over a total range of slightly less than one decade, with respective regression 
coefficients of -0.99 and -0.90. If these slopes were used in accordance with Eq.(|TB|), the 
fractal dimension would he D > 2. Clearly, this cannot be for our planar sets, so that we 
conclude that the procedure suggested by Eq. (|TB|) is not generally valid for the determination 
of fractal dimensions from He scattering data. A further reason is that the intensities 
beyond the specular broadening depend strongly on the hardness of the potential. For 
example, on changing the He/Ag/Pt(lll) potential, a very different dependence, similar 
to that in the inset in Fig.^, is obtained (not shown here). Consequently, one would be 
forced to conclude that the fractal dimension changes when the interaction potential is 
varied. Within a simplistic interpretation, in which D is the mass fractal dimension, this 
is physically unacceptable. If, however. He measures the fractal dimension of the electron 
density contours - see Ref.[|3^ - a dependence on the potential is in fact expected. For 



mass fractals, analytical results have been obtained showing that the scattering intensities 
are self-affine functions, whose scaling exponent is simply related to the dimension of the 
scattering fractal sets.0 
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C. Comparison of Experiments on Ag/Pt(lll) with Scattering Prom Theoretical 

Disorder Models 



In this final section we apply the understanding developed from the analysis of the vari- 
ous disorder models, to interpret He scattering data obtained in the experiments described 
in Sec. [ll A| . To recapture the main details, the system studied in these experiments was a 
submonolayer of Ag at 50% coverage deposited on a Pt(lll) surface at 38K. He scattering 
profiles were taken along the [112] -direction with incident wavenumber of 6.43A~^. The 
experimental data are the circles in Figs-^a-c (the same data is shown in all figures). Our 
purpose in performing the analysis to be described next was threefold: To qualitatively inter- 
pret the available experimental data and thus gain some understanding of the morphology 
of this low-temperature Ag/Pt(lll) phase, to test the insight developed from the analysis 
of the disorder models, and to self-consistently test our empirically extracted potential. 

If one focuses attention on the experimental data in, say, Fig.^, several features (some 
of which observed in the previous sections) stand out: 

1. Sharp specular and Bragg peaks, the latter having "shoulders". 

2. Off-specular structure manifested in oscillations of the intensity. 

3. The absence (or at least masking) of specular peak broadening; instead clear 
maxima at —0.28, 0.37A~^ are present. 

4. An asymmetry between positive and negative AK. 

The last point is an immediate consequence of the experimental scattering geometry, 
which breaks the left-right symmetry (see Sec. [ll A|) . This feature cannot be reproduced in 



the Sudden approximation calculations, in which the He always strikes the surface along its 
normal. The other observations can be used to reach several conclusions on the underlying 
surface structure: 

(0) The large compact clusters model [system (|^)] cannot be ruled out since the broadening 
of the Bragg peaks is one of its main features. 
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(0) The presence of a series of peaks does not agree well with a fractal structure, which 
should exhibit a rather smooth decay of the intensity. Instead, one observes in Fig.pD that 
the scattering pattern from a 50% coverage fractal does contain peaks, for A/^ < 3A~^. 
The reason is that at such high coverage part of the fractal character is lost: adatoms fill 
in the gaps typical of the low coverage fractals of Fig.|^, and the resulting surface structure 
is at first sight not very different from the random one of Fig.Ha. However, the scattering 
intensities are markedly different in Figs-lja and b: the intensity decay in the fractal case is 
far smoother that in the random one, in agreement with the absence of a typical length scale 
(> a) in the fractal. Point ^ then, can be used to eliminate both the fractal and randomly 
adsorbed adatom models: the former due to a lack of peak structure for AK > 3A^^, in the 
latter due to a complete disagreement in peak positions. 

(0) The maxima adjacent to the specular are a feature which was absent in the disorder 
models we studied. They are most likely a result of constructive interference due to a length 
scale of ~ (27r/0.33), i.e. about 4 lattice constants (19/(2. 77v^)) along the [112] direction. 
At present we do not yet have a satisfactory model which contains such a length scale and 
is also in agreement with the other features of the intensity distribution. 

Having essentially ruled out the fractal and random adatoms models, what then is an 
appropriate model which can at least qualitatively explain the features observed in the ex- 
perimental intensities?S The central clue in answering this question is the off-specular peak 
structure. In the absence of agreement with the above models, we attempted a fit with a 
hybrid of the large compact islands and SCC model, namely a narrow distribution of trans- 
lationally random small compact clusters. The advantage of this model is that it turns out 
that the positions and intensities of the off-specular peaks are extremely sensitive to the rel- 
ative numbers and absolute diameters of these clusters. Thus a possible, although definitely 
not final morphological identification, is given in Fig.Hc. The distribution shown is a set of 
hexagonal islands, 7 or 9 Ag atoms in diameter, mixed in a 1:4 ratio at 50% coverage. Upon 
a careful examination it will be noticed that apart from the quantitative location of the 
peaks near the specular discussed above, and the elevated intensity around the first order 
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Bragg, the positions of all other peaks are correctly reproduced (some of these peaks are mere 
shoulders in the experimental data; e.g. at 3.0 and 4.6 A~^). This agreement is not present 
in either the random adatom or fractal cases. On the other hand, the two missing features in 
the intensity pattern of Fig.^ seem likely to be explained by the introduction of some degree 
of positional order in the locations of the SCCs: A separation of about 20-30A between the 
centers of every pair of adjacent clusters would result in the observed peaks, while at the 
same time create a shoulder underneath the Bragg peaks due to a slit-Fraunhofer effect [as 
for system (|^)]. Such an ordered separation is in agreement with the fact that due to the 
size-mismatch between Pt and Ag, a strain is induced by the formation of an Ag island, 
resulting in an effective repulsion between neighboring islands. We observe evidence for this 
in the location of the experimental first order Bragg peak, which is located slightly below the 
position of the corresponding pure Pt(lll) peak, suggesting some degree of lattice ordering 
due to Ag (which has a lattice constant of 2.89A). In addition there are some purely kinetic 
effects determining the size and distribution of islands, namely nucleation and an average 
diffusion length. The generation of a size-dispersed SCCs distribution satisfying the separa- 
tion requirement turns out, however, to be a non-trivial task, and for the present qualitative 
purposes, we limit ourselves to the translationally random model. Beyond the introduction 
of positional order, one should probably also relax the assumption of perfect hexagons: At 
the low surface temperature considered, the mobility of Ag adatoms is extremely small, and 
consequently a more realistic model is probably one in which the hexagons are imperfect at 
least along their perimeter. Since the hexagons are perfect in our model, one observes an 
oscillation in the calculated intensity of Fig.pc, with a period of ~ 0.25A~^; this is likely due 
to interference between adatoms separated by 4-5 lattice constants along the chains making 
up the edges of the hexagons. At any rate, while clearly not being fully quantitative, the fit 
with the narrowly size-distributed SCCs is very promising and bound to be non-accidental. 
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IV. CONCLUDING REMARKS 



The motivation for the study presented in this paper is the rapidly growing interest in 
structures formed on surfaces during epitaxial growth. We demonstrated that various types 
of disordered structures thus formed can be studied quantitatively by He scattering, and that 
a distinction between them can be made on the basis of the different scattering features that 
arise respectively in these cases. This is due to the high sensitivity of the He angular intensity 
spectrum to the surface morphology and electron distribution. In particular, we inverted a 
He/Ag/Pt(lll) potential from experimental cross section measurements, and subsequently 
used this potential to predict from the analysis of experimental angular intensity data the 
presence of a phase of small and compact, narrowly size-dispersed clusters with some degree 
of positional order. The sensitivity to changes in both potential and morphology parameters 
of our fit to the experimental data, as well as support from a more general He-scattering 
studyS and independent STM data,0 provides strong evidence for the existence of this 
phase. Furthermore, it illustrates what we believe to be an important principle in surface 
science: The combination of the mutually complementary experimental (He scattering and 
STM) and theoretical techniques provides an extremely powerful tool for the quantitative 
analysis of complex surface structures. 

Our theoretical findings shed light on several other points: 

1. There are clear and specific differences between the intensity distributions re- 
sulting from the scattering from different classes of disordered structures. Care- 
ful analysis of the relevant features, such as broadening of the specular peak, 
the slope {dlogl /dAK), off-specular Fraunhofer maxima, position and width of 
Bragg peaks and more, can lead a long way to the identification of the type of 
disordered phase present on the surface. 

2. The structure of the complex intensity spectrum obtained by scattering from 
translationally disordered, small compact clusters, can be understood as a sum of 
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the contributions due the scattering spectrum of a single such cluster, and due to 
the disorder. The former is mainly responsible for characteristic Praunhofer and 
rainbow peaks in the large A/T regime, whereas the latter affects the slope of the 
intensity function for small /S.K values and is responsible for the introduction of 
noise and peak "smearing". Conversely, this information can be used to extract 
useful information about the individual cluster size and the cluster ensemble 
statistics. 

3. The claim that for fractals the angular intensity dependence on parallel momen- 
tum transfer follows a universal power law, is not supported by our calculations. 
Instead there is a strong dependence on the parameters of the He/adsorbate po- 
tential. Thus it does not appear to be possible to calculate the fractal dimension 
from a simple scaling analysis of the intensity distribution. 

Recent experimental studies have addressed the different kinds of disordered structures 
formed during 3D epitaxial growth. Specifically, the attenuation of the He scattering spec- 
ular peak during growth was measured, allowing for the distinction between 2D vs. 3D 
growth mechanisms. The results presented here indicate that much more information on 
the structure of disordered multiple ad-layers can be obtained by He scattering, and we are 
currently pursuing such work. 
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APPENDIX 

We derive Eg. (|T5|) . This requires first the evaluation of Eq.(|T^). The task at hand is the 
calculation of the probabilities PT{zj] Ri, R2), to observe a given height difference between 
points located at -Ri and R2 on the surface. Denoting the underlying surface as and the 
atom tops as 1, there exist four possibilities: 

(1) Ri,R2 e with Z{Ri,R2) = 0; 

(2) ^1 e 0, i?2 G 1 with Z{Ri, R2) = h- 

(3) ^1 e 1, ^2 e with Z{Ri, R2) = -h and 

(4) Ri,R2 e 1 with Z{Ri,R2) = 0. 

— * — * 

Denote S'^^ (with e^^ = 0, 1 according as to whether i?j G or i?j G 1), as the probability 
of observing the point Ri in phase e^^ and the point R2 in phase e^^. Then clearly: 

Pr(0;i?i,i?2) = ^oo + ^ii FT{h; Ri, R2) = Soi Fr{-h; Ri, R2) = S,o (17) 

A general formalism for the evaluation of quantities of the type S^^ e~ - e- is given by 
Torquato and Stell (TS).Ell Their result for the case of interest to us is: 

k=l i=l 



N ( \k k 

Sn{Ru .., i?„) = 1 + E ^ j-'-j^^ (^') n 



dxj (18) 



i-na-'^d^.-^.D) 

i=l 

Su{Ri,R2) = 1 - (SiiRi) + S,{R2)) + Soo{Ri,R2) (19) 
5'oi(-Ri, -R2) = 5'i(-Ri) — 5'oo(-Ri, -R2) (20) 

S'io(-Ri, -R2) = 5'i(-R2) — 5'oo(-Ri, -R2) (21) 

where Sn is shorthand for n points in phase (so e.g. Si above means one point in phase 
0), p = N/A is the adsorbate density, x'^ = {xi,...,Xk), and gk is the reduced distribution 
function of k particles, defined as: 

'"(-") = (iV^^'(«') (23) 

Pk (f') = J J ^(^1' ■ ■ ■ dxN- (24) 
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In the case of an isotropic system g2 {xi,X2) = g{r) (where r = \xi — X2I) is the well-known 
radial distribution function. The function m{\Ri — Xj\) is the "particle indicator function", 
defined as: 













if 








if 


1^- 





m{\R-x,\) = { ^ (25) 



Collecting Eqs.(0),(|I§,(|I3), we obtain: 

I{AK) = l^JdRi dR2 e^^^-(^i-^2) X 
(5oo(i?i,i?2) + 5n(i?i,i?2)) + e''"'^5oi(i?i,i?2) + e-''"^5io(i?i,i?2)] ■ (26) 

The last equation, together with Eqs.(p!^)-p^), establish a formal connection between the 
reduced distribution functions of interest in characterization of surface statistics and the 
measurable scattering intensity. 

Eq.(p6|) can next be applied to the scattering from translationally random adatoms. It 
is shown by TS that in this case, since Qk = ^ for all k, one obtains: 

Sn{Ru:,Rn) = e-''''-^^"''\ (27) 

where Vn is the union area of n circles with diameter d, with centers at {Ri, .., Rn}- For 
n = 1, 2 there exist explicit formulae: 

Vi = -nd^ (28) 
4 

V2{R/d) = lnd'- \d^ ||cos-i ^ - f /l^j ^(1 - R/d) (29) 

Here R = \Ri — R2\ and H{x) is the Heavyside step function. Note that in the scattering 
context, d = d{kz), d being the diameter of the cross-section, which depends on the incident 
He wavenumber k^. Using this in Eqs.(^^,(|TU|)-(^T]), and recalling that the present system 
is translationally invariant, we obtain: 

Soo{R/d) = e-P^'^^^l'^') (30) 
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(31) 



Soi{R/d) = Sio{R/d) = - e 



(32) 



What remains is to insert the last expressions into the general formula for the scattering 
intensity, Eg. (PUf) , and go through some algebra and integrations. Defining: 



(3 = pd^ 



7 = 2kz h 
d^ 



(33) 



and: 



ci(/3, 7, A) = ^ (l + 2e-^^(cos(7) - 1)) 



(34) 



C2(/3, 7, e) = 47re(l - cos(7)) e 



(35) 



we obtain: 



/(Ai?) = -i-ci fdRe^^^-^ + ^n, I dRe^'^^-^e-"^^^^"^ 



(36) 



The first integral evaluates to a delta-function, whereas the second simplifies to: 



2ne 2 



Let us now define: 



Jo{RAK)e 



RdR+ Jo{RAK)RdR 

J d 



Jo 



(37) 



then finally the scattering intensity becomes: 



I{AK) = S{AK)c,{P, 7, A) + C2(/3, 7, e) 



g{AKd,P) + 



J_i{AKd) 
AKd 



(38) 
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FIGURES 

FIG. 1. A comparison between experimental and theoretical cross-sections for He scattered from 
an isolated Ag atom adsorbed on a Pt(lll) surface. The experimental values were determined from 
adsorption curves (Eq.(^)) recorded at 38K and extrapolated to — > 0. A fit (solid line) of 
the theoretical cross-sections was used to derive the potential used in this work. 

FIG. 2. Calculated angular intensity distribution along the Tk direction for He scattered from 
an ordered monolayer of Ag on the Pt(lll) surface, with the Ag/Pt(lll) lattice mismatch ignored. 
The orientation of the scattering plane is along the close packed [110] direction in real space. 
The reciprocal space lattice is rotated by 30° so that the Bragg peaks appear at non-close packed 
multiples of An/a, a = 2.77A being the Pt(lll) lattice constant. For the indexing of the Bragg 
peaks, we chose real space unit vectors oi = {a/2,\/3a/2), 0.2 = (—a/2, -v/3a/2), corresponding to 
reciprocal space unit vectors hi = (27r/a)(l, l/\/3), ^2 = (27r/a)(l, — 1/\/3). 

FIG. 3. Dashed line: Calculated angular intensity distribution along the F^^ direction of He 
scattered from a single Ag adatom. Solid line: Same for He scattered from 15% coverage transla- 
tionally random Ag atoms on a Pt(lll) lattice. 

FIG. 4. Dashed line: Calculated angular intensity distribution along the F;^ direction of He 
scattered from a single Ag heptamer. Solid line: Same for He scattered from randomly distributed 
Ag heptamers on Pt(lll). Inset: Top view of a small compact cluster (heptamer). 

FIG. 5. Calculated angular intensity distribution of He scattered from a Pt(lll) surface with 
15% adsorbed Ag atoms with off-lattice translational randomness. Thin line: Sudden approxima- 



tion results. Thick lines: calculation using Eq.(15). Right Inset: Same, using a harder He/Ag 



potential. Left inset: The disorder term g compared to the single adatom term, from Eq.(|15| 
(linear scale). 
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FIG. 6. Calculated angular intensity distribution of He scattered along the Tk direction from 
randomly distributed large compact Ag islands on Pt(lll). Inset: Compact Ag islands produced 
by a KMC simulation, at 15% coverage on a surface of 100 x 100 unit cells and T=500K. The surface 
is a parallelogram with an angle of 60°. Periodic boundary conditions are employed, so straight 
edges usually imply that the cluster is continued at the opposite boundary. 

FIG. 7. Calculated angular intensity distribution of He scattered along the Tk direction from a 
Pt(lll) surface with 15% adsorbed Ag atoms, forming fractal islands. Left inset: KMC simulations 
as in Fig.^, but for fractal islands (T=200K). Right Inset: Power-law regression of the intensity 
(see text for details). 

FIG. 8. A comparison between experimental and theoretical angular intensities for He scattered 
from a Pt(lll) surface with 50% adsorbed Ag atoms in different disorder classes: (a) Translation- 
ally random adatoms, (b) Fractal islands, (c) Size-dispersed compact islands. Insets: Typical 
configurations in each of the disorder classes. The orientation of the scattering plane is along 
the non-close packed [112] direction in real space, corresponding to the close packed Tm direc- 
tion in reciprocal space. Thus Bragg peaks appear at multiples of 47r/(\/3fl) where a = 2.77 A is 
the Pt(lll) lattice constant. Experimental surface temperature: 38K, He beam energy: 21.6meV 
[k^ = QA2>A"^). Solid lines: theoretical intensities. Circles: experimental data. Note that only the 
elastically scattered intensity is shown. The inelastic part was separated in the experiment using 
time-of-flight spectroscopy. 
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